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Abstract. We study quasipolynomials enumerating of proper colorings, nowhere- 
zero tensions, and nowhere-zero flows in an arbitrary CW-complex X, gener- 
alizing the chromatic, tension and flow polynomials of a graph. Our colorings, 
tensions and flows may be either modular (with values in Z/kZ for some k) or 
integral (with values in { — k + 1, . . . , k — 1}). We obtain deletion-contraction 
recurrences and closed formulas for the chromatic, tension and flow quasipoly- 
nomials, assuming certain unimodularity conditions. We use geometric meth- 
ods, specifically Ehrhart theory and inside-out polytopes, we obtain reciprocity 
formulas for the numbers of acyclic and totally cyclic orientations of X. 



This article is about generalizing the enumeration of colorings, flows, cuts and 
tensions from graphs to cell complexes. We begin with a review of colorings of 
graphs. 

Let G — (V, E) be a graph and let A; be a positive integer. A proper k-coloring 
of G is a function / : V — > K, where K is a "palette" of size k and f(v) ^ f(w) 
whenever vw is an edge of G. It is well known that the number Xc{k) °f proper 
fc-colorings is a polynomial in k (the chromatic polynomial of G). Remarkably, 
the numbers xg(~ k) have combinatorial interpretations as well, as discovered by 
Stanley the best known of these is that |xc ( — 1) | is the number of acyclic 
orientations of G. This phenomenon is an instance of combinatorial reciprocity and 
is closely related to, e.g., Ehrhart's enumeration of lattice points in polytopes [19] 
and Zaslavsky's theorems on counting regions of hyperplane arrangements |37j . 

Fix an orientation on the edges of G and let d be the signed incidence matrix of 
G; that is, the rows and columns of d are indexed by vertices and edges respectively, 
and the (u, e) entry is 



If we regard a coloring / as a row vector c = [f{v)] vt =v, then properness says 
precisely that cd is nowhere zero, since for each edge vw, the corresponding entry 
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1. Introduction 
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if vertex v is the head of edge e 
if vertex v is the tail of edge e, 
if v is not an endpoint of e. 
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of cd is f(v) — f(w). Here the palette may be regarded either as the integers 
1, . . . , k or as the elements of TLjkTL = Z^. The number of fc-colorings may then 
be computed by linear algebra and inclusion-exclusion, yielding Whitney's formula 
for the chromatic polynomial (see, e.g., |36[ §5.3]). 

The cut and flow spaces of G are respectively the row space and kernel of its 
boundary matrix <9, regarded as a map of modules over a ring A (typically either Z, 
R, or Zfc). Cuts and flows arise in algebraic graph theory and are connected to the 
critical group and the chip-firing game: see, e.g., pQ or [3TJ Chapter 14]. Over Z, 
we consider in addition the space of tensions, or integer vectors of whom some 
multiple is in the cut space. (Over a field or a finite ring, all cuts are tensions.) For 
a generalization of cuts and flows to cell complexes, see [15] . 

If we regard the graph G as a one-dimensional simplicial complex, the matrix d is 
just the boundary map from 1-chains to 0-chains. Accordingly, we can replace the 
graph by an arbitrary d-dimensional CW-complex X and define cellular colorings, 
flows, tensions and cuts in terms of its top cellular boundary map d : Cd(A;Z) — > 
C d -i(X;Z). 

In Section [3j we study the functions 

X*x(k) — number of proper Z^-colorings of X, 
T x (k) = number of nowhere-zero Zfc-tensions of X, 
(p* x (k) = number of nowhere-zero Z^-flows of X. 

When X is a graph, these are all polynomials in k (for the tension and flow polyno- 
mials, see [35]); in fact, they are specializations of the Tutte polynomial of X. For 
arbitrary cell complexes, we show that the modular counting functions are always 
quasipolynomials, and find sufficient conditions on X for them to be genuine poly- 
nomials (building on the work of the first author and Y. Kemper [2]). We describe 
two avenues toward such results, using colorings as an example (the other arguments 
are similar). First, if d contains a unit entry, then pivoting the matrix d there cor- 
responds to the topological operations of deleting a facet or deformation-retracting 
it onto a neighboring ridge, and gives rise to a deletion-contraction recurrence for 
modular colorings. Provided that this "shrinking" process can be iterated, it follows 
by induction that Xx(^) * s a polynomial. The second approach is a linear-algebraic 
inclusion-exclusion argument, which produces a closed quasipolynomial formula for 
X*x(k) (Theorem 3.4 1, which is easily seen to be polynomial if every column-selected 
submatrix of d has no nontrivial invariant factors (a weaker condition than total 
unimodularity) . 

In Section [4j we study the integral counting functions 

Xx (k) = number of proper A"-colorings of X, 

i~x(k) = number of nowhere-zero A-tensions of X, 

fx(k) = number of nowhere-zero A-flows of X 

where A is the palette {— k + 1, — k + 2, . . . , k — 1} C Z. When X is a graph, these 
are all polynomials in k as was shown by Kochol [55[ [53] . These can be regarded 
as Ehrhart functions, enumerating lattice points in inside-out polytopes [1]. While 
we cannot in general find explicit combinatorial formulas for the integral counting 
functions, we can interpret them combinatorially in terms of acyclic and totally 
cyclic orientations of X (generalizing the graph-theoretic definitions), and obtain 
reciprocity theorems: for instance, |xa(~1)| is the number of acyclic orientations 
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of X. Using these results, we obtain reciprocity theorems for the modular counting 
functions in Section [5] by interpreting them as Ehrhart quasipolynomials of certain 
modular analogues of inside-out polytopes. 

The authors thank Hailong Dao, Yvonne Kemper and Jacob White for numerous 
helpful suggestions. 

2. Background 

The symbol N will denote the nonnegative integers. For n G N, we set [n] := 
{1, . . . , n}, and for a < b G Z, we set [a, b] :— {a, a + 1, . . . , b — 1, e}. 

We start by defining colorings, tensions and flows in the general context of CW- 
complexes. We assume familiarity with the basic machinery of CW-complexes, 
including simplicial and cellular homology. Unless otherwise stated, we will follow 
the notation and terminology of Hatcher [33] . 

Throughout the paper, X will denote a cell complex of dimension d. We assume 
that X is pure, i.e., every maximal cell has the same dimension. (We lose no gen- 
erality by this assumption; it is analogous to assuming that graphs are connected.) 
The symbols F = {fi, . . . , f m } and R = {ri, . . . , r„} will always denote the sets of 
d-cells and (d — l)-cells of X; these are called facets and ridges, respectively. We 
allow the possibility that F = (even though, strictly speaking, the dimension of 
such a complex is less than d). For an Abelian group A, the chain groups Cd(X; A) 
and Cd~i(X;A) can be identified with the free modules A F and A R respectively. 
The boundary map Ci{X) — > Cj_i(X) is denoted 3a,i or d Ai , and the coboundary 
map Ci-i(X) —> Ci(X) is d* A i . We simplify this notation by dropping the super- 
script when the complex X is clear from context; moreover, if A or i (or both) are 
not specified, the reader should assume by default that A = Z and i = d. The map 
d is represented by an n x m matrix with columns indexed by F and rows by R, 
and entries d r f given by the degrees of the cellular attaching maps. The (reduced) 
homology groups are 

H i {X;A)=kerdlJlmd Ati+1 . 

The cellular matroid of X is the matroid M(X) represented by the columns of 
dR y d- (For a general reference on matroids, see, e.g., [29].) In the case that X is a 
graph, the cellular matroid is just the usual graphic matroid; more generally, a set 
J C F is independent in M(X) if Hd{Xj; Z) = (where Xj is the CW-subcomplex 
of X generated by J U R) . It is often convenient to use matroid language when 
speaking of facets or sets of facets of X. For instance, a loop of X is a facet whose 
boundary is zero, and a coloop is a facet whose boundary is linearly independent 
of the boundaries of the other facets. 

Definition 2.1. An A-coloring of I is a vector c = (ci,...,c„) € A R . The 
coloring c is proper if cd contains no zero entries. An integral k-coloring of X is a 
Z-coloring c such that c, G [— k + 1, k — 1] for every i. 

Definition 2.2. The A-flow space of X is the A-module 

Flow A (A) := kcrcU U f = H d (X;A). 

An integral k-flow is a Z-fiow (u>i, . . . , w m ) such that Wi G [— k+ 1, k— 1] for every i. 
In the graph setting, an enumeration theory of fc-flows was initiated by Kochol [2S] ; 
see also [SJ H2] . 
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Definition 2.3. An orientation of X is a sign vector e £ {±1} F that labels each 
facet of X with a positive or negative sign. 

A flow w £ Flowz(A) is e-nonnegative if £iWi > for every facet fi £ X. The 
terms e-positive, s -negative, etc., are defined analogously. We write ew to denote 
the vector obtained by componentwise multiplication. A cell complex X with given 
orientation s is acyclic if it has no nontrivial e-nonnegative flows. It is totally cyclic 
if for every facet /j, there exists a flow w that is e-nonnegative and SiWi > 0. An 
equivalent condition is that X has a strictly e-positive flow (since if . . . , u/ m ) 
are e-nonnegative flows with SiW^ > for each i, then ^\e«;W is a strictly e- 
positive flow). 

Let A F be equipped with the standard inner product (a/)/gir • (6/)/eF = 
^2feF a fbf- For every Z-submodule £> C A F , define B- 1 = {a e A F : a • 6 = 
for every &€_£?}. 

Definition 2.4. The A-cut space and A-tension space are 



These are submodules of A F ; note that Ten A {X) = (Cut A (X)- L )- L D Cut^(X) and 
that Tenz(X)/ Cutz(A) is isomorphic to the torsion summand of H d (X;Z). An 
integral k-tension of A" is a Z-tension -0 such that tpi £ [—k + 1, k — 1] for every i. 

When A is a graph, the numbers of nowhere-zero Z^-cuts, flows and tensions are 
polynomial functions of k [35] and the polynomials for cuts and tensions coincide. 
In the cell complex setting, we will see that they are quasipolynomials. A quasipoly- 
nomial is a function q(k) of the form q(k) = c n (k) k n + c n _i(fc) k n ~ x + ■ ■ • + Co(fc), 
where Co(fc), . . . , c n (k) are periodic functions of k; the least common period of 
co(fc), . . . , c n {k) is the period of (Thus a polynomial is a quasipolynomial 

of period 1.) 

Proposition 2.5. Let A be a field or a finite commutative ring. Then Cut^(A) = 



Proof. We will show that (B ) = B, where B is a submodule of a finitely gener- 
ated free A- module M — A r . This is elementary if A is a field. If A is finite, then 
the inner product identifies M r with its Pontryagin dual M* = Hom(Af, Q/Z). 
(Usually the Pontryagin dual is defined as Hom(M, C x ), but we wish to write 
the group law additively.) For any A-submodule (equivalently, Abelian subgroup) 
B C M, there is a map M* — > B* , given by restriction, that is surjective (since 
any homomorphism B* — > Q/Z can be extended to a homomorphism M* — > Q/Z) 
and has kernel B 1 - . We therefore have an isomorphism M*/B* = B- 1 ; there- 
fore |M* = |M|/|B| = B^- (where the first equality follows from Pontryagin 
duality). In particular | (i? -1 -)- 1 - 1 = \B\, and since B C (B- 1 )- 1 we must have equal- 
ity. Taking B — Cutfc(A), we obtain the desired equality of the cut and tension 
spaces. □ 

The following useful fact is immediate from the foregoing definitions. 

Proposition 2.6. The map c <— > ip :— cd is a surjective 1-linear map from color- 
ings to cuts. The coloring c is "proper if and only if the cut ip is nowhere zero. 



Cut A (X) 
Ten j4 (A) 



Rowspace^ <9, 
Flow A (A)- L . 



Ten A (X). 
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All of our definitions are natural extensions of the usual notions of cuts, flows 
and colorings for graphs. A flow of a graph is an element of its cycle space: an 
assignment of weights to the oriented edges so that the total flow into each vertex 
equals the total flow out. A coloring is a function on the vertices. For graphs, cuts 
and tensions are synonymous: the minimal supports of nonzero elements of the cut 
space are its bonds (the minimal edge sets whose deletion increases the number of 
components, or equivalently the cocircuits of the graphic matroid). The definitions 
of the cut and flow space of a cell complex are the same as those of [5] and [TS] 
(with different notation). 

2.1. Unimodularity conditions. Let M £ z nxm be a matrix of rank p. We say 
that M is totally unimodular (TU) if and only if the determinant of each square 
submatrix is —1, 0, or 1. Totally unimodular matrices arise in several contexts, 
such as integer linear programming [311 Chapter 8] and as representations of reg- 
ular matroids [29, Chapter 6]. A totally unimodular cell complex is one whose top 
boundary map is TU. In topological terms, X is TU if and only if for all subcom- 
plexes Z C Y C X with dim!" = d and dimZ — d— 1, the relative homology group 
Hd-i{Y, Z; Z) is torsion-free, as proved by Dey, Hirani and Krishnamoorthy [TTI 
Theorem 5.2]. For example, every graph is TU. 

The following facts about TU matrices may be well known, but we have not 
found an explicit reference in the literature, so we include short proofs here. For a 
general reference on TU matrices and their properties, see, e.g., |3U Chapter 19]. 

Lemma 2.7. Let M be a n x m TU matrix and B S 1 n . Then the system of 
equations MX = B has a rational solution in X if and only if it has an integer 
solution. 

Proof. The "if" direction is clear. For the "only if" direction, let Xq be a rational 
solution. Then X is a solution to the system of linear inequalities AX < C, where 



A = 



M 




B 


-M 


c = 


-B 


I 


P 


-I 




_ -Q J 



are matrices in block form, where I is the m x m identity matrix and P, Q are 
obtained by rounding the coordinates of A up and down, respectively, to the nearest 
integers. Then A is a TU matrix [3TJ eqn. (43), §19.4], and the solution space 
of {X : AX < C} is a nonempty polytope contained in the affine linear space 
{X: MX ~ B}. By [3"T1 Thm. 19.1], the vertices of P are all integral, so any such 
vertex is the desired solution. □ 

Lemma 2.8. Let m > n and let M be a n x m TU matrix. Then there exists a 
(m — n) x m TU matrix Z such that the rows of Z are a It-basis for ker A. 

Proof. Write M = [B\C]. Assume, without loss of generality, that B is a column 
basis for M, so that it is also a Z-basis for Colspace z (M). Note that B € SL n (Z). 
By applying repeated pivot operations, we can transform the first r columns into 
the identity matrix; this is equivalent to multiplying M on the left by B , so 
the resulting matrix is [J|B _1 C]. By [521 Lemma 2.2.20], this matrix is TU. Let 
Z = [{B~ 1 C) T \ - I] e Z( m ~ n ) xm . Then the rows of Z are a vector-space basis 
for kerM, and moreover Z is TU because every minor of Z equals some minor of 
[IIB^C]. □ 
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In order to state our enumerative results for as general a class of cell complexes 
as possible, we will consider various weakenings of the TU condition, as we now 
explain. 

Definition 2.9. The matrix M is quasi-unimodular (QU) if the greatest common 
divisor of its maximal-rank minors is 1. It is strongly quasi-unimodular (SQU) if 
every column-selected submatrix of M is QU. 

We recall some standard facts about finitely generated Abelian groups. There 
exist unique positive integers ax, . . . ,a p , the invariant factors of M, such that 
ai\a 2 \ ■ ■ ■ | o p and cokerM := Z n / Colspace(M) = Z n -^©0f =1 Z ai . These numbers 
can be computed as ai — gi/gi-i, where gi denotes the gcd of the i x i minors of M 
(putting g = 1). In particular the invariant factors of a matrix and its transpose 
are identical. 

Define a quasipolynomial function of M by 

p 

(1) 1 (M,k) = IJgcdCM*). 

1=1 

Knowing the function 7(M, k) is equivalent to knowing the isomorphism type of 
the group coker M. Moreover, if we denote by M ® the matrix M with entries 
regarded as elements of Z^, then 

(2) |ker(M®Z fe )| = k n - p j(M,k). 

If X is a cell complex of dimension d, we put j(X, k) = j(d* , k), where d* — d'^(X) 
is its top coboundary matrix. Then we have 

The product a^a-i ■ ■ ■ a{ equals the greatest common divisor of the i x i mi- 
nors of M. Thus the boundary matrix of a CW-complex X is QU if and only if 
Hd-i(X;Z) is torsion-free, and it is SQU if and only if Hd-i(Y;Z) is torsion-free 
for every subcomplex Y C X generated by a subset of facets. 

2.2. Deletion and contraction. Let d be a n x m matrix over A, and let r G [n] 

and / € [m] . Define a matrix d' by deleting the f th column from d. If d r f is a unit 
in A, then we may in addition define a matrix d" by pivoting at d r f. 

• For each i =^ r, replace the row di by di — d r . 

• Delete the row and column corresponding to r and / respectively. 
That is, the entries of d" are 

(4) % S =0iJ-% L 8ri 

If d is the (integer) cellular boundary matrix of a cell complex X and d r t = ±1, 
then these operations carry topological meaning: d' and d" are respectively the 
boundary matrices of the deletion X' — X \ /, and the "contraction" X/rf — 
X/(r, /) obtained by deformation-retracting / onto r. For instance, if X is a graph 
then these complexes are the usual deletion and contraction of the edge /. More 
generally, the column matroids of & and d" are the deletion and contraction of 
that of d HH Proposition 3.2.6]. 
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Total unimodularity of an integer matrix is preserved by pivot operations. There- 
fore, any TU complex can be repeatedly contracted, choosing the facet and ridge 
arbitrarily at every step, to eventually yield a cell complex whose boundary map is 
either zero or empty. More generally: 

Definition 2.10. Let X be a cell complex and let (r, /) be a ridge-facet pair. Then 
X is shrinkable via (r, /) if d r f = ±1. It is iteratively shrinkable (ISH) if d is zero 
or empty, or if it is shrinkable at some pair (r, /) and X \ f and X/rf are both 
iteratively shrinkable. 

Thus every TU complex is ISH; in fact, any non-loop facet in an ISH complex 
can be chosen as the first one to shrink. Moreover, every ISH complex is QU, 
because if successively shrinking via the pairs (r l5 fi) . . . , (r s , /„) produces a zero 
or empty matrix, then the submatrix with rows r\ , . . . , r s and columns /i, . . . , f s 
has determinant ±1. 

To summarize, we have the implications 

TU =>> SQU =► QU 

and 

TU =*> ISH => QU. 

On the other hand, the conditions SQU and ISH are incomparable. For instance, 
any matrix invertible over Z is SQU but need not be ISH, while a matrix with a 
full-rank identity submatrix is ISH but need not be SQU. 

2.3. Inside-out polytopes. A famous theorem of Stanley [32] gives a combina- 
torial interpretation for the values of the chromatic polynomial Xg(^) a ^ negative 
integers k; in particular, |xg( — 1)1 is the number of acyclic orientations of G. Simi- 
lar interpretations have been found for the values of the integral and modular flow 
and tension polynomials of graphs [4[ [5J [9] . In Sect ions [4] and [5j we will extend these 
results to cell complexes using lattice-point geometry, specifically Ehrhart theory 
(see, e.g., [5]) and inside-out polytopes [1], which we review here. 
Let S be a subset of K™. Its Ehrhart function is 

(5) E s (k) := \Z n nkS\ 

where kS denotes the k th dilate of S for k € 1i>\. If Q is a rational convex polytope 
(i.e., the convex hull of finitely many points in M. n with rational coordinates), then 
its Ehrhart function EqIU) has the following properties [3J Theorem 3.23, Corollary 
3.15]: 

• Eq is a quasipolynomial in k of degree dimQ. In particular, this allows us 
to evaluate Eq at non-positive integers, even though |5| is only defined for 
positive integers. 

• The period of Eq divides the denominator of Q (i.e., the smallest k such 
that k~ lr L n contains every vertex of Q). 

• E Q (0) = 1. 

We are going to be interested in the Ehrhart functions of various Boolean combi- 
nations of rational convex polytopes. Note that for any such Boolean combination 
5 1 , the above results automatically imply that Es{k) is a quasipolynomial and can 
thus be evaluated at non-positive integers. Of special importance are the Ehrhart 
quasipolynomials of open polytopes: If Q° is the relative interior of a rational convex 
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polytope Q, then Eqo (fc) is a quasipolynomial and the famous Ehrhart-Macdonald 
reciprocity theorem [31 Theorem 4.1] states that 

(6) (-lf™QE Q (-k) = E Q o(k). 

Let W — {Hi, . . . , H m } be a central hyperplane arrangement in E". The corre- 
sponding inside- out polytope [¥] is 

m 

C:= (-1,1)"\|J^- 

j'=l 

Let . . . , i?° be the open regions (i.e., connected components) of C, and for each i, 
let Ri denote the topological closure of R° . Note that each Ri is a convex polytope 
and that C is the disjoint union of the R°. Therefore, the Ehrhart quasipolynomial 
of C is given by 

s 

E c (k) :=^^o(fc). 

i=l 

We also define the closed Ehrhart quasipolynomial of C as 

s 

E c (k) :=^E Ri (k). 

i=l 

For any positive integer k, Ec(k) counts lattice points z G (Z™ n fcC) with a 
multiplicity given by the number of regions Ri that z is contained in. Note 
that Ec(0) = 5Zi = i -Ejii (0) = s. Our motivation for defining the closed Ehrhart 
quasipolynomial of C is ([6| which gives us 

(-l) n E c (-k) = E c {k). 

Motivated by this reciprocity, we will sometimes call the Ehrhart quasipolynomial 
of C the open Ehrhart quasipolynomial of C . 

3. Modular Counting Functions 

Throughout the rest of the paper, X will denote a cell complex of dimension d, 
with facet and ridge sets F,R of cardinalities m,n respectively, and top cellular 
boundary matrix d of rank p. 

In this section, we study the enumeration of Z^-colorings, tensions and flows on a 
cell complex X. The modular chromatic, tension and flow functions are quasipoly- 
nomials in k. Subject to certain unimodularity assumptions, they satisfy deletion- 
contraction recurrences and Whitney-style inclusion-exclusion formulas, and can 
be recovered (essentially) from the Tutte polynomial of the cellular matroid of X. 
Many of the results generalize well-known results about the chromatic, tension, and 
flow polynomials of a graph; see, e.g., [351 §5-3]. 

3.1. Modular colorings. Recall that an A-coloring of X is a row vector c = 
[c r ]re-R such that cd is nowhere zero. More generally, a coloring of a matrix M is 
a row vector c such that cM is nowhere zero. 

Definition 3.1. The modular chromatic function of X is 

Xx{k) '■= number of proper Z^-colorings of X. 
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When X is a graph, this is just its chromatic polynomial. We will show that in 
general, the modular chromatic function is a quasipolynomial, and it is a polynomial 
when X is SQU. 

If X has no facets, then x*x( k ) = ^™ (because all colorings are proper), while 
if d has at least one zero column, then Xx{k) ~ 0. X is a graph, the second 
condition says that X has a loop.) Moreover, if P is any n x n matrix that is 
invertible over Z, then replacing d with Pd does not change the modular chromatic 
function of X (even though this operation does not in general make topological 
sense). 

Theorem 3.2. Let r and f be a ridge and facet of X , respectively, such that 
d rS = ±1, and let X' = X\f and X" = X/rf. Then 

X*x( k ) = X*x'(k)~X*x"( k )- 

Proof. Every proper coloring c of X is a proper coloring of X\f (because if cd = 0, 
then cd'j = cdj for all j 6 F \ /). Conversely, if d — (cj)j 6 .R is a proper coloring 
of X' and ip — cd, then rpj ^ for all j ^ f , and c is a proper coloring of X if and 
only if ipf =/= 0. 

Suppose that ipf = 0. Then we claim that this proper coloring of X' restricts to 
a proper coloring c" — (c\, . . . , c r , ■ ■ ■ , c n ) of X" (where the hat denotes removal). 
Indeed, for j ^ /, we have by Q 



(c" • d") 3 = £ fa - |^) 



i£R\r iGR\r 




C r d rJ C-rdrf ~ ^ Cjdj/ 

ie.R / 



J rf 



= to - ihf o. 

On the other hand, we claim that every proper coloring c" = (c±, . . . , c r , . . . , c„) 
of X" can be extended uniquely to a proper coloring d of X' with cdf =0. Indeed, 
there is a unique value of c r for which c<9/ = 0, namely 

(V = — — Cidif , 

rf i£R\r 

and then, for all j ^ f, we have (again using (|4 

i£R ieR\r ^ r ' ' Ti i^r ieR\r 

which is nonzero by assumption on c". Accordingly, the correspondence between d 
and c" is a bijection. □ 

Remark 3.3. This argument does not work for integral fc-colorings. That is, if 
d r f = 1 and c G Z R ^ r ^ is a proper fc-coloring of X/rf with —A: < c, < fc for every 
i, then we may define c r = X^e-R-fr} c i&ifi but there is no reason to expect that 

— k < Cf < k. 

We now consider the case that / is a coloop of the cellular matroid. Again, 
suppose that d r f = 1 and that / and r are represented by the first column and first 
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row of d. Pivoting at the (1, 1) position (a series of Z-invertible row operations) 
gives a matrix of the form 



" 1 


V 





d" 



in which v is Q-linearly dependent on the rows of d" . (Note that this condition 
is always satisfied when X is SQU .) If v is in fact a "L-linear combination of the 
rows of d" , then we can perform further Z-invertible row operations to obtain the 
matrix 



" 1 








ffl 



whose colorings are precisely the row vectors of the form [c\, . . . , c r ] such that c\ £ 
Zzk \ {0} and [02, . . . , c r ] is a coloring of X" . Therefore, under these assumptions, 
we have 

(7) X*x(k) = (k-l) X *x»(k). 

We now consider how to express the modular chromatic function as a Whitney- 
style sum over subsets of facets, or as an evaluation of the Tutte polynomial. For 
a set J C F of facets, recall that Xj is the CW-subcomplex of X generated by 
J U R. Let dj be its boundary matrix, which is the submatrix of d consisting of 
the columns indexed by J. Let p(J) = rank<9j; note that p is the rank function of 
the cellular matroid A4(X). 

Theorem 3.4. We have 

X*x(k) - ^(-l)l J lfc"-l J l|iJ rf pO;Z fc )| - ^(-l^lfc^M^fc) 

JCF JCF 

where 7 is as defined in ([T]). In particular, the modular coloring function x*x(k) 
is a quasipolynomial in k whose period divides the 1cm of the set of all invariant 
factors of column- selected submatrices of d. 

Proof. By inclusion-exclusion and ^ we have 

X*x(k)= ^(-l) |J| |{ceZf : C a, 7 = 0}| 

JCF 

= |ker Zfc d}\ 

JCF 

= £(-l)l J lfc»-l J l|H d (Xj;Z fc )| 

JCF 

= J2(-V lJl k n - r l(dj,k), 

JCF 

proving the first assertion, from which the second follows as a special case. The 
third assertion follows from the formula for the Tutte polynomial as a corank-nullity 
generating function [TU]. □ 

Corollary 3.5. Suppose that X has no zero columns and d is SQU (for instance, 
if X is TU). Then 

Xx(k) = ^(-l)l J lfc"-"( J ) = (-l)"k n -"T x (l-k,Q) 

JCF 
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where T is the Tutte polynomial of M.(X). In particular, Xx(^) * s a m onic poly- 
nomial in k of degree n that is divisible by k n ~ p , whose coefficients alternate in 
sign, and whose k 12 ^ 1 -coefficient is the number of equivalence classes of facets of X 
under the relation that two facets are equivalent if the corresponding columns of d 
are scalar multiples of each other over Q. 

Proof. These assertions follow from the observation that Xx(^) i s J us t the charac- 
teristic polynomial of the hyperplane arrangement whose normal vectors are the 
columns of d; see, e.g., [321 §3.11.2]. □ 

Example 3.6. Let X be the real projective plane M.P 2 with its standard cell 
structure (one cell in each of dimensions 0, 1, 2). The top boundary map is d — [2], 
so a modular K-coloring is an assignment of a color c € 7L k to the 1-cell such that 
2c ^ mod k. Therefore, the modular coloring function is 

a. , , \ k — 1 if k is odd 
X*x(k) = { = fc-gcd(fc,2). 

Ik — 2 it fc is even 

Example 3.7. Let X = Ajv.,2 be the d-dimensional skeleton of a simplex on N 
vertices, so that m — (Jlj) and n — (^). This cell complex is QU (because X 
is homotopy-equivalent to a wedge of d-dimensional spheres), but is not TU in 
general; for example, Ae,2 contains subcomplexes homeomorphic to M.P 2 . When 
d = 1, A„,,i is a complete graph, whose chromatic polynomial is well known to be 
k(k — l)(fc — 2) • • • (k — n + 1). A general formula for Xx(^) m terms of N and d is 
not apparent. 

3.2. Modular cuts and tensions. 

Definition 3.8. For k £ N, let T^-(k) be the number of nowhere-zero Zfc-tensions. 
The function t x : N — > Z is called the modular tension function of X. 

Proposition 3.9. We have 

T *(k) = -= X *(k) = 4 TTX*x(k). 

XK ' \H d (X;Z k )\ XXK ' j(X,k) Xxy 1 
Proof. Recall that the cut space and the tension space are identical over a finite 



base ring (Proposition 2.5). Therefore, the nowhere-zero tensions are precisely the 
vectors of the form c d for some proper coloring c, and each such tension arises in 
this way from exactly | ker d* \ colorings. □ 



By combining Proposition |3.9| with the results of Section |3.1| we obtain the 
following ways to compute r x {k). 



Corollary 3.10. (1) (Inclusion-exclusion formula) By Theorem 3.4 we have 



rUk) = 1 V (-l) ]Jl k m -W\H d (Xj:Z k )\ 



l{X,k) 



(-i)i ^kp-p^^x -j, k). 
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(2) (Deletion-contraction) Suppose that d r f = I. Let X' = X — f and X" = 
X/rf. Then 

r* x ( k) = ig^zoi 

\H*(X;Z k )\ X \H*(X;Z k )\ 
^X',k) l{X",k) 
j(X,k) kTx ' {k) 7 (X,fc) Tx " [k) - 

If in addition f is a coloop and the row vector [d r j]j£F\{f}_is a Hi-linear 
combination of the vectors \dij\j£F\{f} fori € i?\{r} ; then fcfh implies that 

r x (k) = (k-l)r x „(k). 

(3) (Tutte-polynomial specialization) If X is SQU and T is the Tutte polynomial 
ofM(X), then 

r x (k) = (-l) r T x (l-k,0). 

3.3. Modular flows. 

Definition 3.11. For k € N, let <p x (k) denote the number of nowhere-zero Z^-flows 
on X . The function ip x : N — > Z is called the modular flow function of X. 

Theorem 3.12. Suppose that X is shrinkable via the ridge-facet pair {r,f). Let 
X' = X\f and X" = X/rf. Then 

tp* x (k) = tp* x „(k) - (Px'(k)- 

Moreover, if w is a flow on X, then w\x/ r f * s a fl ow on X/rf and if w is a flow 
on X/rf then there exists a flow w on X with w\x/ r f = w ' ■ 

Proof. Without loss of generality, assume that d r f = 1. By definition, a nowhere- 
zero Zfc-flow on the contraction X" is a vector (w%, . . . , w} , . . . , w m ), with all G 
Zfc \ {0}, such that 

uiidi = mod df\ 
that is, there is a unique scalar Wf such that 

w f d f = Widi ■ 

If uif ^ then (w±, . . . ,w m ) is a nowhere-zero flow on X, while if Wf = then 
(u>i, . . . ,w), . . . , w m ) is a nowhere-zero flow on X' . 

Meanwhile, every nowhere-zero flow w on X or on X' gives rise to a nowhere-zero 
flow on X" (in the former case, by omitting wA, □ 

If X has no facets then f* x {k) = 0, and if X has a loop / then 

tp* x (k) = kip* XXf (k) . 

Thus by induction on the number of facets, we have: 

Corollary 3.13. If X is ISH, then (p* x (k) is a polynomial in k. 

Theorem 3.14. We have 

<p*xW = ^i-ir-^k^-^hixj,^. 

JCF 
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In particular, the modular flow function (p* x (k) is a quasipolynomial in k whose pe- 
riod divides the 1cm of the set of all invariant factors of column-selected submatrices 
ofd. 

If X has no zero columns and d is SQU (for instance, if X is TV), then 

JCF 

In particular, in this case the modular flow function is a monic polynomial of de- 
gree m, and 

<p* x (k) = (-l) m -"T x (0,l-k) 
where T x denotes the Tutte polynomial of M(X). 

Proof. Fix fceN. By inclusion-exclusion and ([2]), we have 

= E (- 1 ) m_|J 'l{^ G ke rz t d: w f = for all / € F\ J}\ 



JCF 

E 

JCF 

E 

JCF 



(-ir-| J l|ker Zfc dj\ 
(-l)™-\ J \k\ J \-p( J ^(Xj,k) 



as desired. The third assertion follows from the formula for the Tutte polynomial 
as a corank-nullity generating function [10] . □ 

Remark 3.15. The columns of d represent an arithmetic matroid in the sense of 
d'Adderio and Moci [Til §3.2]. The corresponding arithmetic Tutte polynomial is 

M x {x,y)=Y J \T(H d (Xj;Z))\(x - ly-^Hy - 1)^-^ 

JCF 

where T denotes the torsion summand. In [15], the same authors specialize the 
arithmetic Tutte polynomial to enumerate colorings and flows on certain generalized 
edge-weighted graphs. Our modular chromatic, tension and flow functions do not 
appear to be recoverable from M x , because the quasipolynomial j(Xj, k) specifies 
the actual group structure of T(H d (X; Z)) rather than just its cardinality. On the 
other hand, the modular chromatic and tension functions can be obtained from the 
Tutte-Grothendieck class of the arithmetic matroid represented by d [5UJ §7]- 

4. Integral Counting Functions and Reciprocity 

In this section, we extend various reciprocity results for integral counting poly- 
nomials defined on graphs to the more general integral counting quasipolynomials 
defined on cell complexes. To this end, we use the general framework of inside-out 
polytopes as developed by Beck and Zaslavsky [H[5]. In particular, the reciprocity 
theorem for the integral chromatic quasipolynomial of a cell complex, Theorem |4. 5 [ 
generalizes Stanley's reciprocity theorem for the chromatic polynomial of a graph 
|32) . In this context, Proposition 4.4 generalizes Greene's geometric insight about 
acyclic orientations of graphs [55 . The reciprocity theorem for the integral ten- 
sion quasipolynomial, Theorem 4.8 generalizes reciprocity results by Chen [llj 
and Dall [16] . Finally, the reciprocity theorem for the integral flow quasipolyno- 
mial, Theorem |4.11[ generalizes the integral flow reciprocity theorem by Beck and 
Zaslavsky [S]. 
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4.1. Integral colorings and acyclic orientations. Recall that an integral k- 
coloring of X is an element of (— k, k) R (note that there are 2k — 1 available colors 
for an integral £;-coloring) . A coloring is proper if c d is nowhere zero. 

Definition 4.1. The integral coloring junction of X is 

Xx(2fe — 1) := number of proper integral fc-colorings of X. 

If we regard a coloring c as a lattice point in Z R , then the proper /c-colorings are 
precisely the lattice points in kC, where C is the inside-out polytope (— 1, 1) R \ / Hx , 
where T-Lx is the arrangement of hyperplanes {Hf : f € F} normal to the columns 
of d. Therefore, the integral coloring function satisfies 

(8) X x(2k - 1) = E c (k) = (-l) n E c (-k) . 

Definition 4.2. Recall that an orientation of X is a sign vector e € {—1, 1} F . We 
will abuse notation by also writing e for the m x m diagonal matrix with entries 
Sf for / G F. Reorienting X by e means replacing d with cte, i.e., multiplying each 
column of d by the corresponding entry of s. 

Also recall that an orientation e is acyclic if it admits no nontrivial nonnegative 
flows; that is, every nonzero element of Flow(A, e) = ker(cte) has at least one 
positive and at least one negative entry. An equivalent condition is that des =/= 
for every s € N F \ {0}. 

Let c be an integral coloring and let e be an orientation. We say that the 
pair (c, e) is compatible if the vector cde has strictly positive entries, i.e., Ef = 
1 whenever J2reR c r^rf > 0, and Ef = — 1 whenever J2r£R c r^rf < 0. Every 
coloring has 2 a compatible orientations, where a = \{f € F: J2 r eR c rdrf = 0}|. 
In particular, a coloring is proper if and only if it has exactly one compatible 
orientation. 

Proposition 4.3. Let c be a proper coloring of X and let ip = cd be the corre- 



sponding nowhere-zero tension (see Proposition 2.6). Define an orientation e by 
Then e is acyclic. 

Proof. It suffices to show that des ^ for every s E N F \ {0}. Indeed, for such s, 

cdes = ip ■ (s f Sf) f< =F = ^f £ f s f + 
feF 

because ipfEf > ioi every /, and it follows that des ^ 0. □ 

Proposition 4.4. There is a bijection between regions ofHx and acyclic orienta- 
tions of X . 

Proof. Every open region Q° of M. R \ U/gF Hf ^ s unbounded, hence contains some 
lattice point, which is a proper coloring and therefore gives rise to an acyclic orien- 
tation by the previous lemma. Moreover, all colorings in Q° give rise to the same 
acyclic orientation, denoted £q, in this way. If two regions are distinct, then they 
lie on opposite sides of at least one hyperplane Hf, and so the corresponding orien- 
tations will differ in (at least) their f th position. Accordingly, it remains to show 
that every acyclic orientation e of X has a compatible integral coloring. This is a 
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direct application of Gordan's theorem from linear programming |13[ Exercise A. 3, 
p. 334] to the definitions of acyclic orientation and compatibility. □ 

Theorem 4.5. Let X be an arbitrary cell complex. The number of pairs (e, c), 
where e is an acyclic orientation of X and c is a k-coloring compatible with e, 
equals {— l) n Xx(— 2fc — 1)- In particular, the number of acyclic orientations of X 
is (-irxx(-l). 

Proof. Equation Q may be rewritten as 

E c (k) = {~l) n E c (-k) = (-l)" X x(-2fc - 1) 

where C is the inside-out polytope (—1, 1)\(J Hx- Meanwhile, Ec(k) is the number 
of pairs (R, c) , where c is a fc-coloring and R is a closed region of Hx such that 
c £ R. Furthermore, for each region R, the lattice points in R° are precisely the 
colorings compatible with the orientation corresponding to R, which implies the first 
assertion of the theorem. The second assertion follows by setting k = 0, because 
the only 1-coloring of X is the zero vector, which lies in every closed region. □ 

4.2. Integral tensions and acyclic orientations. Recall that the tension space 
Ten(X) is the orthogonal complement of the flow space Flow(X) = ker<9. In 
particular, the tension space is a saturated Z-submodule of Z F (i.e., if some nonzero 
scalar multiple of -0 is a tension, then so is ip), which is equivalent to the statement 
that Ten^(A) is a summand of Z F . 

Recall also that an integral tension (resp., a fc-tension) is a tension ip £ 1 F (resp., 
(-k,k) F ). 

Definition 4.6. The integral tension function of X is 

Tx(2fc — 1) := number of nowhere-zero /c-tensions of X. 

Let B be the Boolean arrangement of coordinate hyperplanes in M. F . Define the 
inside-out polytope 

T:= ((-l,l) F nTen(X))\B. 
Then the nowhere-zero /c-tensions are precisely the lattice points in kT. 

If X has a loop /, then Ten(X) is contained in the corresponding coordinate 
hyperplane. In this case, there are no nowhere-zero tensions; the inside-out polytope 
T is empty; rx(2fc — 1) = 0; and X has no acyclic orientations. We assume 



henceforth that X has no loops (but note that Proposition 4.7 and Theorem 4.8 
below still hold trivially). 
We have 

dim(T) = dim(Ten(X)) = m - dim(Flow(X)) = p. 
Accordingly, the closed and open Ehrhart quasipolynomials of T satisfy 
(9) T X (2k-l) = E T (k) = {-lYE T (-k) 

where the second equality comes from Ehrhart-Macdonald reciprocity. 

A tension tp and an orientation e are called compatible if, for all facets /, we have 
e fi ) f > 0. That is, £/ = 1 whenever tjjf > 0, and £/ = — 1 whenever -0/ < 0. Note 
that every tension i/j has 2 a compatible orientations, where a = |{/: ip* = 0}|; in 
particular, a tension is nowhere zero if and only if it has exactly one compatible 
orientation. 

Proposition 4.7. There is a bijection between regions of Ten(A) \ B and acyclic 
orientations of X. 
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Proof. Every open region R° of Ten(X) \ B is unbounded, hence it contains some 
lattice point, which is a nowhere-zero tension and therefore gives rise to an acyclic 
orientation by Proposition |4.3| Moreover, all tensions in R° give rise to the same 
acyclic orientation sr in this way. If two regions are distinct, then they lie on 
opposite sides of at least one coordinate hyperplane, and so the corresponding 
orientations will differ in that coordinate. Accordingly, it remains to show that 
every acyclic orientation e of X has a compatible nowhere-zero tension. By Propo- 
sition |4.4| every acyclic orientation e has a compatible proper coloring c. Then 
by the definition of compatibility tp — cdx is a nowhere-zero tension compatible 
with e. □ 

Theorem 4.8. Let X be an arbitrary cell complex. The number of pairs (s,if)), 
where e is an acyclic orientation of X and ?p is a k-tension compatible with e, 
equals (— l) p Tx(— 2fc — 1). In particular, the number of acyclic orientations of X is 

(-i)M-i). 



The argument is the same as that in the proof of Theorem 4.5 
4.3. Integral flows and totally cyclic orientations. 
Definition 4.9. The integral flow function of X is 

(px{2k — f) := number of nowhere-zero fc-flows of X. 

As before, let B be the Boolean arrangement of coordinate hyperplanes in M. F . 
Define the inside-out polytope 

W := ((-1, i f n Flow(A)) \ B . 

Then the nowhere-zero £;-flows are precisely the lattice points in kW. 

If &x has a column that is linearly independent from all others (i.e., the matroid 
represented by X has a coloop), then Flow(X) is contained in the corresponding 
coordinate hyperplane. In this case W = and ipx(2k — 1) = 0, and X has no 



totally cyclic orientations (and Theorem 4.11 below is trivially true). Henceforth, 



we assume that X has no coloops, in which case 

dim(VF) = dim(Flow(X)) = m — p. 

Accordingly, the closed and open Ehrhart quasipolynomials of W satisfy 

(10) fxC2k - 1) - E w (k) = (~l) m - p E w (-k) 

where the second equality comes from Ehrhart-Macdonald reciprocity. 

Every lattice point in W gives rise to a totally cyclic orientation (see Defini- 
tion 2.2 1, namely the vector of signs of its coordinates. Moreover, lattice points in 
the same region of W give rise to the same totally cyclic orientation. On the other 
hand, for any totally cyclic orientation e, there is a strictly e-positive flow that is 
a lattice point in W, and all e-positive flows belong to the same region of W . In 
summary we obtain a proposition analogous to Proposition |4.7| 



Proposition 4.10. There is a bisection between regions ofF\aw(X)\B and totally 
cyclic orientations of X . 

We say that a fc-flow w and a totally cyclic orientation e of a cell complex X are 
compatible if ew > 0. 
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Theorem 4.11. Let X be an arbitrary cell complex. The number of pairs (e,w), 
where e is a totally cyclic orientation of X and w is a k-flow compatible with e, 
equals (— l) m ~ fl tpx{— 2fc — 1). In particular, the number of totally cyclic orientations 
ofXis(-l) m -P<p x (-l). 

Again, the argument is the same as that in the proof of Theorem |4.5| 

5. Modular Reciprocity 

In this section, we generalize various reciprocity theorems for modular counting 
polynomials defined in terms of graphs to the case of modular counting quasipoly- 
nomials defined in terms of cell complexes. To this end, we follow the geometric 
approach developed by Breuer and Sanyal .9!, by realizing the counting functions 
as Ehrhart quasipolynomials of a disjoint union of polytopes. This approach makes 
important use of the correspondence between regions and orientations that we de- 
veloped for the integral case in the previous section. In particular, the modular flow 
reciprocity theorems, Theorem |5.2| and Corollary |5.5[ generalize the modular flow 
reciprocity theorem for graphs by Breuer and Sanyal [5j . See Chen and Stanley [T2] 
for a related reciprocity theorem. The reciprocity theorems for modular tensions, 



Theorem 5.6 and Corollary 5.8 generalize reciprocity results by Chen [llj and 



Breuer and Sanyal [5]. The reciprocity theorem for modular colorings, Theorem 5.9 
generalizes the classic reciprocity for the chromatic polynomial by Stanley [35], see 
also Beck and Zaslavsky [3]. 

In order to obtain reciprocity theorems for the modular flow and tension func- 
tions, we again interpret them as lattice points in inside-out polytopes, following 
the approach in [9]. 

If v = (vi, . . . , Vk) is a vector, we write 

supp(u) — {i: Vi 0}, zero(w) = {i: Vi = 0}. 

5.1. Modular flow reciprocity. Recall that the modular flow function <f* x (k) 
counts the nowhere-zero Z^-flows of X. For an integer vector b £ Z R , define 

P d {b) := {w £ [0, 1] F : dw = b} , P$(b) := {w G (0, 1) F : dw = b) . 

We call b feasible if Pg(b) ^ 0. Let B = B(d) denote the set of all feasible vectors, 
and define 

Pa := [J Pa(b), P° d := |J P° d {b)- 

beB beB 

Then the modular flow function of X is the Ehrhart function of Pq , that is, 

<p* x (k) = Ep S (k) = \z F nkP$\. 

Note that for all b € B the dimension of Pg(b) is m — p, since the intersection 
of a translate of Flow^(X) with the open set (0, 1) F has the same dimension as 
Flow^(X). Therefore, we can apply Macdonald-Ehrhart reciprocity to obtain 

(-1)— />25 P o(-fc) = £(-i)— P23 p . (6) (-fc) 

beB 
beB 

= E Pa (k) 
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which gives us 

(11) (-i) m -"ip* x (-k) = \z F n kP d \ . 

Note that if C is an n x n invertible integer matrix, then 

B{C^d) = {b£Z R : {w E (0, l) F : C~ l dw = b} ^ 0} 
= {b£Z R :{we (0, 1) F : dw = Cb} ^ 0} 
= C- 1 {a e Z R : {w e (0,1) F : dw = a} ^ 0} 
= C^Bid) 

and therefore 

(12) P c - i9 = [J P c -x a (6)= □ P c - ia (C- 1 a)= □ P„(a) = P 8 , 

beB(C~ 1 d) a£B(d) a£B(d) 

so we will henceforth drop the map d from the notation. 



Example 5.1. Let d = 



1 2 

2 4 



. Then 



P(0,0) = {(0,0)}, P(l,2) = £((1,0), (0, 1)), 

P(2,4) = £((1, |), (0, 1)), P(3, 6) = {(1, 1)}, 

where £(p, q) denotes the closed line segment from p to q. The feasible vectors are 
B{d) = {(1,2), (2,4)}. The sets Pg(b) are the open line segments shown in blue in 
the leftmost figure below. The number of interior lattice points in the k th dilate is 
the quasipolynomial 

k — 1 if A; is odd, 
k — 2 if k is even. 




Figure 1. (a) The "modular flow inside-out polytope" of Exam- 
ple 5.1 (b), (c) Its fc-dilates for k — 3 and k = 4. The black dots 



are the points of Pg, counted by ip* x (k). 



Recall from Section |4.3| that X is coloop-free if and only if it has at least one 
totally cyclic orientation. In addition, X has no coloops if and only if ker9 is 
not contained in any coordinate hypcrplane, i.e., it contains a nowhere-zero vector. 
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Indeed, if d has an nowhere-zero nullvector, then each column of d can be solved for 
as a linear combination of the other columns. Conversely, if d has no coloops, then 
for each facet / there is a nullvector with a nonzero element in the f th position, and 
taking a generic real linear combination of these nullvectors gives a nowhere-zero 
vector in FIowr(A). 

Our modular analogue of Theorem |4.11| can be most easily expressed through 

!w is a Z/j-flow on X and 
(w, a): o : zero(u)) -> {— 1, 1} extends 
to a totally cyclic orientation of X 

Theorem 5.2. Suppose that X has no coloops. Then 

{-xr--<p* x {-k) = \* k {x)\. 

Proof. For w £ Z F n kPg, define h(w) — (w, a), where w — w mod k and 

a f = 

We will now show that h is a bijection. By (llll) this will complete the proof. 




First, to show that h is well-defined, we observe on the one hand that, by con- 
struction, w is a Zfe-flow on X because dw — kb' for some b' £ 1 R . On the other 
hand, a totally cyclic orientation e of X with a — e| zcr o(iD) can be found as follows. 
Let b £ 7L R such that w £ kP(b). By construction 

kP°(b) = (FIowk(X) + w) n (0, k) F . 

In particular, kP°(b) is a relatively open, non-empty set of dimension to — p. Since 
X has no coloops, the flow space is not contained in any coordinate hyperplane 
(that is, any hyperplane in the Boolean arrangement B), and the intersection of any 
coordinate hyperplane with FIowr(A) has dimension to — p — 1. By the identity 

kP°{b) \{B + w) = ((Flow R (X) \ B) + w) n (0, k) F , 

it follows that kP°{b) \(B + w) can be obtained from kP°(b) by removing finitely 
many sets of dimension m—p—1. This implies in particular that kP°(b) \ (B + w) is 
non-empty. Let w £ kP°(b)\(B+w) and v := w— w. As both w, w £ F1owr(X)+u>, 
it follows that v £ FIowr(X). As w g (B + w), it follows that v £" B. Thus, 
v is a nowhere-zero R-fiow on X. Let e denote the sign vector of v. By the 



reciprocity theorem for integral flows (Theorem 4.11) we know that e is a totally 
cyclic orientation of X. Moreover, for / £ zero(w)), either Wf = 0, when u)f > wt, 
or else Wf = k, when wt < Wf, and in either case Sf = o~f. 

It is immediate from the definition that h is injective. To show surjectivity, 
suppose that w is a Z^-flow on X and a is a function zero(w) —¥ { — 1, 1} such that 
there exists a totally cyclic orientation e of X with a — e| ZG ro(is)- Define w £ Z F 
by 

if / zero(w), 
if / £ zero(w) and <jf = — 1, 
if / G zero(u') and crj = 1. 

where lift (a;) denotes the unique integer in {1, 2, . . . , k — 1} that reduces to cc mod- 
ulo k. Then cku = kb for some & € 1 R (because w is a Z^-flow) and w £ Z F n [0, k] 




F 



So all that remains to show is that Pg(b) ^ 0. By Theorem 4.11 there exists a 
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nowhere-zero real now v with sgn(u) = e; we may choose v to have arbitrarily small 
magnitude. We claim that w = w — v G P° } (b). Indeed, first, we have 

dw = dw — 8v = kb — = kb . 

Second, if / ^ zero(w) then < Wt < k. If Wf = then sgn(wj) = —1 and Wf > 0, 
while if Wf = k then sgn(uj) = 1 and Wf < k. This proves the claim and completes 
the proof that h is surjective. □ 

Corollary 5.3. If X has no coloops, then \ip* x (—l)\ is the number of totally cyclic 
orientations of X . 



Proof. If we set k = 1 in the definition (13), we see that the vector w must be the 
zero vector, and a — e is a totally cyclic orientation of X. □ 

Remark 5.4. Suppose that X has a coloop f\. Let C be an invertible r x r 
matrix with an initial set of columns that form a Z-basis for the saturation of 
^{/2j ■ • • ; fm}- Then C~ 1 d has the form 



ai 






8 


O-p-l 




a p 











Note that d is the boundary map of the complex X/f obtained by collapsing the 
cell / to a point, and w = (w\, . . . , w m ) is a flow on X if and only if (u>2, ■ ■ • , w m ) 
is a flow on X/f and W\a,i = mod k for all i = 1, 2, . . . , p. The number of values 
of W\ for which the second condition holds is gcd(ai, . . . , a p , k). Therefore, 

f*x{k) =gcd(ai,...,a p ,fc) <p* X / f (k). 

Note that if X is a TU complex, then for any set S of facets the contraction X/S 
is well defined. In this case the modular flow reciprocity theorem for complexes can 
be phrased in exactly the same way as the modular flow reciprocity theorem for 
graphs. 

Corollary 5.5. Let X be a totally unimodular complex without coloops. Then 
\<p*x(—k)\ counts the number of pairs (Hi, a) where w is a Zk-flow on X and a is a 
totally cyclic orientation on X/supp(w). 



Proof. To derive the corollary from Theorem |5.2[ all we have to show is that for any 
Zfc-flow w and any a : zero(w) — > {—1, 1} the following statements are equivalent: 

(1) a extends to a totally cyclic orientation of X. 

(2) a is a totally cyclic orientation of X/ supp(u)). 

Suppose that supp(u)) = {fx, ■ ■ ■ , ft}- By the observation after Definition 2.10 
we can find a sequence of ridges r\,...,ri such that the contractions Xq = X, 
Xi = X /(rx,fx), Xi = Xi_x/{n, fi), ■ ■ ■ , X t = X/ supp(w) are well-defined. 

(JTJ) implies |2j): Let e be a totally cyclic orientation of X such that e| zcro ( tB ) = cr. 
Then by Proposition 4.10 there exists an £-positive E-flow v on X. We claim that 
w lzcro(to) is a CT-positive flow on X/supp(w). Clearly, v\ zexo tw) is cr-positivc. By 
Theorem 3.12 we know that if w|jrw I) „,/ 4 _ 1 is an M-flow on X^i, then w|_F\/ l! .../ i 
is an M-flow on X t . By induction we obtain that v\ ZCTQ ^ is an M-flow on Xg = 
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X/ supp(«)). Thus, by Proposition 4.10 we conclude that a is a totally cyclic 
orientation on X/ svpp(w). 

([2]) implies ([lj: Let a be a totally cyclic orientation of Xj supp(w). Then there 



exists a er-positive IR-flow vg on Xg = X/ supp(?Z>). By Theorem 3.12 
inductively construct K-flows vi, . . . , vq = v on Xg, -X^_i, . . . , Xq = X such that 
sgn(^| zero ( lB )) = a. (Note that these flows need not be nowhere-zero.) 

We now construct an IR-flow w on X such that zero(u>) = zero(u)). First, inter- 
pret w as a vector in 7L F by lifting all entries to their canonical representatives in 
Z. Then dw — kb for some b £ l Z R , Second, choose a lattice point z £ 1 F such that 

2| ZC ro(ts) = and dz = b. 
(Such a lattice point must exist becaus e this linear system of equations has a rational 



solution, namely -rid, and by Lemma 2.7 it must also have an integral solution.) 
Finally, let w = w — kz, so that dw = 0. Note that wj 7^ whenever / £ supp(u)), 
because by construction, Wf is not a multiple of k. 

Now, v+Xw is an M-flow on X for any A £ R. Because wr = for all / £ zero(w), 
we can choose A such that v + Xw is nowhere zero and has the same sign as v on 



all facets / £ zero(w). By Proposition 4.10 there exists a totally cyclic orientation 



e compatible with v + Xw, and we have e| ze ro(«j) = c as desired. □ 

5.2. Modular tension reciprocity. Recall that the modular tension function 
rjj-(fc) gives the number of nowhere-zero Z^-tensions of X, or equivalently (by 



Proposition 2.5| the number of nowhere-zero Z^-cuts. In this section we again use 
lattice-point methods to obtain a reciprocity theorem for modular tensions. The 
theorems and proofs are similar, and in many places dual, to those of Section [5.1| 
Let S be a set of vectors that forms a basis for ker d as a Z-module. Let [S] be 
the matrix whose rows are the vectors in S. Then the tension space 

Ten R (X) = {V £ R F \ [S]ip = 0} 

is of dimension p. For each b £ Z R , define 

Q s (b) := {ip £ [0, 1] F : [<% = 6}, Q° s (b) := £ (0, 1) F : [S\1> = b}. 

Call b feasible if Qg(b) 7^ 0. Let B(S) denote the set of all feasible vectors, and 
define 

Qs:= U Qs{b), Q° s := [J Q° s (b). 

beB(S) beB(S) 

Then the modular tension function of X is the Ehrhart function of Q5, that is, 

T* x (k)=E Q% (k) = \Z F nkQ° s \. 

For each b £ B, the dimension of Q^(b) is p, since the intersection of a translate of 
Tenj{(X) with the open set (0, 1) F has the same dimension as Ten^X). Therefore 
we can apply Ehrhart-Macdonald reciprocity to obtain 

(-l)^Qi(-fe) = X)(-l)'iSgo(6)(-*) 

b£B 
b£B 

= E Qs {k) 
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which gives us 

(14) (-l) p r* x (-k) = \Z F nkQ s \. 

As before, for ip £ Z F , let zero(?/>) denote the set of all / £ F with ipf = 0. 

The construction of Q is independent of the choice of S, for the following reason. 
If S' is another Z-basis for ker d, then there is an invertible Z-matrix C of size m — p 



such that C[S] = [S'], and a calculation similar to (12) shows that 
Q°cs(b) = Q° S {C- x b% Q° cs (Ca) = Q° s (a), 

B{S) = C^BfjOS), CB{S) = B(CS), 

Qs = Qcs, Q°s = Qcs- 

For this reason, henceforth we write Q instead of Qs . 

If X has no loops — i.e., d has no zero columns — then a generic real linear combi- 
nation of its rows gives a nowhere-zero vector in CutjR(X), hence in Tena(X). Thus, 
if X has no loops, then TenR(cc) is not contained in any coordinate hyperplane. 



To state the modular analogue of Theorem 4.8 we define 



ip is a Zfe-tension on X and 
(15) &k(X) := { (V>, cr) : o : zero(^) {-1, 1} extends 

to an acyclic orientation of X 

Theorem 5.6. Suppose that X has no loops. Then 

(-l)'r£(-fc) = |* fc (J0|. 

Proof. For ip e Z F f] kQ, define h(tp) = (0, cr), where if) — ip mod k and 



°~f 



-1 ifV/ = 0, 
1 if ip f = k. 



We will show that h is a bijection. By (14) this will complete the proof. 



First, to show that h is well-defined, we observe on the one hand that, by con- 
struction, ^ is a Zfc-tension on X because d^i = kb' for some b' £ 1 R . On the other 
hand, an acyclic orientation e of X with a — £| zcro (^) can be constructed as follows. 
Let b £ 1 R such that ip £ kQ(b). By construction 

kQ°{b) = (Ten K (X) + V) H {0,k) F . 

In particular, kQ°(b) is a relatively open, non-empty set of dimension p. Since X 
has no loops, the tension space is not contained in any coordinate hyperplane, and 
the intersection of any coordinate hyperplane with TenR(JT) has dimension p — 1. 
We have 

kQ°{b)\{B + i>) = ((Ten R (X)\£0 + V)n(O,fc) F , 

and so kQ° (b)\(B+ip) can be obtained from kQ°(b) by removing finitely many sets 
of dimension p — 1. This implies in particular that kQ°{b) \(B + ip) is non-empty. 
Let £ kP°(b) \(B + ip) and v:=ip- tp. As both ip, ip £ Ten R (A) + w, it follows 
that v £ Tenu(A). As ip £" (B + ip), it follows that v £" B. Thus, v is a nowhere-zero 
K-tension on X. Let e denote the sign vector of v. By the reciprocity theorem 



for integral tensions (Theorem 4.8) we know that e is an acyclic orientation of X. 
Moreover, for / £ zero( - 0), either ipf = 0, when ipf > ipf, or else ipf — k, when 
ipf < ipf, and in either case Sf = a f. 
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It is immediate from the definition that h is injective. To show surjectivity, 
suppose that ip is a Z^-tension on X and a is a function zero(^) — > { — 1, 1} such 
that there exists a totally cyclic orientation s of X with a = £| zcro (^,)- Define 
ip € Z F by 

riift(^) if /0zero$), 
-0/ := s if / € zero(V>) and 07 = — 1, 

[fc if / € zero(V') and cr/ = 1. 

Then [<S]-0 = /cfe for some 6 € Z R (because ip is a Z^-tension) and ■0 G Z F n [0, fc] F . 



So all that remains to show is that Q°(b) ^ 0. By Theorem 4.8 there exists a 
nowhere-zero real tension v with sgn(v) = e; we may choose v to have arbitrarily 
small magnitude. We claim that ip = ip — v € Q°(b). Indeed, first, we have 
[5]^ = - [S]v = kb - = kb. Second, if / g zero(^) then < -0/ < ft. If 
V>/ = then sgn(i>y) = —1 and -0^ > 0. while if ipt = k then sgn(w^) = 1 and 
ipf < k. This proves the claim and completes the proof that h is surjective. □ 

Corollary 5.7. If X has no loops, then |tJ-(— 1)| is the number of acyclic orien- 
tations of X . 



Proof. If we set k = 1 in the definition (15), we see that the vector ip must be the 



zero vector, and a = e is an acyclic orientation of X . □ 

Corollary 5.8. Let X be a totally unimodular complex without loops. Then fc)| 
counts the number of pairs (w, a) where w is a "L^-tension on X and a is an acyclic 
orientation on X \ supp(w). 



Proof. To derive the corollary from Theorem |5.6| all we have to show is that for 
any Z^-tension w and any a : zexo(w) — ¥ {—1, 1} the following statements are 
equivalent: 

(1) er extends to an acyclic orientation of X. 

(2) cr is an acyclic orientation of X \ supp(u)). 

Suppose that supp(u)) = {fx,..., ft}. Define Xj := X \ {fx, . . . , f} for i = 
0,...,£. Thus I = Xandl £ = X\ supp(w). 
We will need the following two observations: 

(a) Let Y be a cell complex and let / be a facet of Y. If ip is an R-tension on 
Y, then ip\y\f is an R-tension on Y \ f. Conversely, if ip' is an R-tension on 
Y \ f, then there exists an R-tension ip on Y such that ip\y\f = V 1 '- This is 
easy to derive from the fact that ip is an K-tension on Y if and only if there 



exists c € R such that ip = cdy; see also Proposition 2.6 
(b) If X is TU, then [S] can be chosen to be totally unimodular, by Lemma 2 

Now we can turn to the proof proper. 

(JTJ) implies |2|): Let e be an acyclic orientation of X such that e\ zeTO (w) = cr. 



Then by Proposition 4.7 there exists an e-positive ]R-tension ip on X. We claim 
that ^lzero(iD) is a tr-positive tension on X\supp(w)). Indeed, ip\ zelo (w) is cr-positive, 
and by observation Q, if ip restricts to an R-tension on then it restricts to an 
R-tension on Xi. By induction ip\ ze m(w) is an R-tension on Xi> = X\supp(u)). Thus, 
by Proposition 4.7 we conclude that a is an acyclic orientation onI\ supp(zD). 
([2]) implies (TFT Let a be an acyclic orientation of X \ supp(iD). Then by 4.7 



there exists a tr-positive R-tension ip e on X e = X\supp(w). By observation Q, we 
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can inductively construct K-tensions ipiy i^i-ii "00 — ip on %l> ■ ■ ■ > Xq — X such that 
sgn( - i | zero ( lD )) = a. Note that these tensions need not be nowhere-zero, although 
ip is nonzero on zero(u>). 

Next, we construct an R-tension w on X such that Wf = if and only if / G 
zero(w). Such a w can be obtained as follows. First, we interpret w as a vector in 
1 F by lifting all entries to their canonical representatives in Z. Then [S]w — kb for 
some b G 7L R . Now we choose a lattice point z G "L F such that z| ZC ro(tD) = and 
[S]z = b. Such a lattice point exists, because jriZS is a rational solution to this linear 
system. Since [5] is totally unimodular by observation Q , this system also has an 
integral solution by Lemma |2.7| t Finally, we let w = w — kz whence dw = 0. Note 
that Wf ^ whenever / G supp(u)) because in this case Wf is not a multiple of k. 

Now, ip + Xw is an K-tension on X for any A G R. Because w/ — for all 
/ G zcro(iD), we can choose A such that tp + Xw is nowhere zero and has the same 
sign as tp on all facets / G zero(u)). Let e be the acyclic orientation compatible 
with if] + Xw, which exists by 4.7 We then have e\ zem {w) — a as desired. □ 



5.3. Modular coloring reciprocity. In analogy to the definition of the set ^>k{X) 



in (15), let 



c is a Z/j-coloring on X and 
(16) Ck(X) := { (c, er) : a : zcro(c<9) — > {—1, 1} extends 

to an acyclic orientation of X 

Theorem 5.9. Suppose X has no loops. Then 

(-l) n X * x (-k) = \C k (X)\. 

Moreover, \x*x(~ 1)1 equals the number of acyclic orientations of X. 

Proof. Consider the map / : Ck{X) -> ^k(X) defined by /(c, a) = (cd, a). This 
map is surjective because every tension is a cut over (Proposition |2.5[ ) and the 
cut space is by definition the rowspace of d. By pL for each ip G Ten^(X) we have 

ir^V^HIkerzja^Hfc'^^fc). 

Therefore, 

\c k {x)\ = k n -p 1 {d 1 k)\^k{x)\ 



(-l) p fc"- p 7(<9, k)r* x {-k) (by Thm. 5.6 1 



(-l) n Xx (-k) (by Thm. 3.91 



The second assertion now follows from Corollary 5.7 □ 



6. Further directions 

We conclude with some possibilities for future research. 

(1) We do not know whether there exists a non-SQU cell complex X whose 
modular chromatic function x*x(k) is an honest polynomial (that is, a 
quasipolynomial of period 1). This would be an instance of "period col- 
lapse"; see, e.g., [35]. For this to occur, all non-monomial terms in Theo- 
rem [32] would have to cancel out. 

(2) For a graph G, the modular flow polynomial V'gW counts the number of 
flows not merely over Z^, but over any Abelian group A of cardinality k. 
This is easily seen to be false for general cell complexes; the matrix d — 
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[2] has two flows over Z4 but four over Z2 © Z2. What if we impose a 
unimodularity condition on X? 

(3) Kook, Reiner and Stanton [57] gave a formula for the Tutte polynomial of 
a matroid as a convolution of tension and flow polynomials. Breuer and 
Sanyal [9j §5] used the Kook-Reiner-Stanton formula, together with reci- 
procity results to give a general combinatorial interpretation of the values 
of the Tutte polynomial of a graph G at positive integers; see also [3D] and 
[SJ Theorem 3.11.7]. Do these results generalize to cell complexes whose 
tension and flow functions are polynomials? 

(4) In the graph case, the geometric setup has proven extremely useful for es- 
tablishing bounds on the coefficients of the chromatic polynomial (cf. Hersh 
and Swartz |24j ) and the tension and flow polynomials, in particular in the 
modular case (cf. Breuer and Dall [5]). Moreover these geometric construc- 
tions are closely related to Steingn'msson's coloring complex [33] [7J. Can 
these methods be extended to the case of counting quasipolynomials de- 
fined in terms of cell complexes? In particular, what are good bounds on 
the coefficients of the quasipolynomials introduced in this paper? 
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